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Temperature-dependent gap equations in the SU(4) model of high-Tc superconductivity are de- 
rived and analytical solutions are obtained. Based on these solutions, a generic gap diagram de- 
scribing the features of energy gaps as functions of doping P is presented and a phase diagram 
illustrating the phase structure as a function of temperature T and doping P is sketched. A special 
doping point P q occurs naturally in the solutions that separates two phases at temperature T = 0: 
a pure superconducting phase on one side (P > P q ) and a phase with superconductivity strongly 
suppressed by antiferromagnetism on the other (P < P q ). We interpret P q as a quantum phase tran- 
sition point. Moreover, the pairing gap is found to have two solutions for P < P q : a small gap that is 
associated with competition between superconductivity and antiferromagnetism and is responsible 
for the ground state superconductivity, and a large gap without antiferromagnetic suppression that 
corresponds to a collective excited state. A pseudogap appears in the solutions that terminates at 
P q and originates from the competition between d-wave superconductivity and antiferromagnetism. 
Nevertheless, this conclusion does not contradict the preformed pair picture conceptually if the 
preformed pairs are generally defined as any pairs formed before pairing condensation. 

PACS numbers: 71.10.-w, 74.20.Mn, 74.25.Dw, 74.72.-h 



I. INTRODUCTION 

The properties of of high-temperature (high-Tc ) super- 
conductivity in copper oxide materials pose a strong chal- 
lenge to the theoretical understanding of superconductiv- 
ity in strongly correlated many-body systems. The main 
question concerning the high-Tc superconductivity phase 
diagram is the transition between the antiferromagnetic 
(AF) and superconducting (SC) phases, which is domi- 
nated by anomalous properties commonly attributed to 
a pseudogap in the spectrum . Partly because of the 
complicated material properties and the nature of the ex- 
periments performed, different experiments seem to em- 
phasize different aspects of the system, leading to differ- 
ent and sometimes contradictory possibilities for explana- 
tion of observables ■ Many models have been proposed 
in the last two decades but there is no consensus as to 
the origin of the pseudogap properties. Still, one cannot 
answer the question: what is the real phase diagram in 
cuprates? 

Describing collective motion in a strongly correlated 
many-body system in terms of single particle degrees of 
freedom is not always feasible. In the cuprate systems 
exhibiting high-T c superconducting properties, there is a 
substantial point of view (see, for example, Ref. 0) that 
the many-body correlations are so strong that the dy- 
namics can no longer be described meaningfully in terms 
of individual fundamental particles. However, collective 



'Electronic address: ysun@nd.edu 
^Electronic address: guidry@utk.edu 



* Electronic address: clwu@phys.cts.nthu.edu.tw 



motions in such quantum many-body systems are often 
governed by only a few collective degrees of freedom. 
Once these degrees of freedom are identified and prop- 
erly incorporated into a model, calculations may become 
feasible and, more importantly, the physics may become 
transparent. One systematic approach to this program 
of identifying the relevant collective degrees of freedom 
in a many-body system is the method of dynamical sym- 
metries. 

Examples of using the dynamical symmetry approach 
to study many-body problems may be found in nuclear 
physics 0, 0, molecular physics 0, 0, and particle 
physics Recently, we have developed an SU(4) dy- 
namical symmetry model aimed at understand- 
ing the mechanism that leads to high-T c superconductiv- 
ity in cuprates. Theoretical models employing algebras 
and groups have been used in condensed matter physics 
(for recent examples, see Refs. However, to our 
knowledge, the powerful dynamical symmetry methods 
that we employ here have not been applied systemati- 
cally in this field. 

In this paper we derive and solve the temperature- 
dependent gap equations in the SU(4) model. We demon- 
strate that this model can provide important insight 
into the puzzling issues associated with cuprate super- 
conductors. In the next section we outline the SU(4) 
model (a more detailed description can be found in Refs. 
IH1 El)- Temperature-dependent gap equations of 
the SU(4) model are derived in section III. Analytical 
solutions for the gap equations at T = are obtained 
in section IV, and a generic gap diagram (energy gaps 
versus doping at T = 0) predicted by the SU(4) model 
is presented. The solutions of the gap equations at finite 
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temperature are given in section V and these are used to 
construct a generic phase diagram in section VI. Finally, 
a summary is given in section VII. Technical details con- 
cerning construction, derivation, and solution of the gap 
equations are given in three appendixes. 



II. THE SU(4) MODEL 

It is by now widely agreed that in cuprates, the mech- 
anism responsible for superconductivity is closely related 
to the unusual antiferromagnetic insulator properties of 
their normal states. There are also compelling arguments 
that the pair mechanism leading to high-T c superconduc- 
tivity does not correspond to ordinary BCS s-wave pair- 
ing. Phase-sensitive experiments indicate that the SC 
phase of most cuprates has <i-wave-like pairing symme- 
try (at least for the hole-doped compounds), and this 
is also supported by photoemission experiments, which 
show the existence of nodal points in the quasiparticle 
gap. 

Such observations argue strongly for a theory based 
on continuous symmetries of the dynamical system that 
is capable of describing more sophisticated pairing than 
found in the simple BCS picture, and capable of unify- 
ing SC and AF collective modes and the corresponding 
phases. Then such fundamentally different physics as SC 
order and AF order can emerge from the same effective 
Hamiltonian as global variables (e.g., doping and tem- 
perature) are varied. 



A. Basic ingredients of the SU(4) model 

The basic assumption of the SU(4) model is that the 
configuration space is built from coherent pairs formed 
from two electrons (or holes) centered on adjacent lattice 
sites. [The SU(4) model is particle-hole symmetric [l4^. 
For cuprate superconductors one is generally interested 
in hole-doped compounds but more general applications 
of SU(4) can deal with either electrons or holes. Here- 
after, unless specified explicitly, we use "electrons" to 
reference either electrons or holes.] In cuprates the co- 
herent pairs are believed to exhibit d-wave orbital sym- 
metry 15], and we assume a coexistence of two kinds 
of coherent pairs in a minimal model: the spin-singlet 
(D) and the spin-triplet (ir) pairs. We adopt the d-wave 
pairs with structure defined in Refs. [13,113 as the basic 
dynamical building blocks of the SU(4) model 
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where ct i creates an electron of momentum k and spin 
projection i,j = 1 or 2 (=f or |), 

g(k) = (cos k x — sinfcy) 



is the (i-wave form factor, and Q = (ir, it, it) is an AF or- 
dering vector. These pair operators, when supplemented 
with operators of particle-hole type Qij and SV,- , 
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constitute a 16-element operator set that is closed under 
a U(4) D U(l) x SU(4) algebra if the condition 

g(k) ps sgn (cos k x — cos k y ) 

is imposed. In Eq. (0), fl is the maximum number of 
doped electrons that can form coherent pairs, assuming 
the normal state (at half filling) to be the vacuum. The 
U(l) factor in U(l) x SU(4) is associated with charge- 
density waves and is independent of the SU(4) algebra. 
The charge-density waves can be excluded in the symme- 
try limit and in the following discussion we shall restrict 
attention to the SU(4) subgroup. The 15 SU(4) genera- 
tors are related to more physical operators through the 
linear combinations 
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where S is the spin operator, Q the staggered magnetiza- 
tion, 7?t (7?) the vector form of the creation (annihilation) 
operator of spin-triplet pairs, and n the electron number 
operator. It will also sometimes prove useful to replace 
the number operator n with the charge operator M, de- 
fined through 

M= ~(Sn+S 22 ) = |(n-«). 

It has been demonstrated that this SU(4) algebra 
defines the minimal symmetry implementing a unified 
description of antifcrromagnetism and d-wave supercon- 
ductivity that is consistent with Mott-insulator proper- 
ties for normal states. 



B. The group structure and model Hamiltonian 

The SU(4) group has three dynamical symmetry group 
chains: 

D SO(4) x U(l) D SU(2) s x U(l) 
SU(4) D SO(5) D SU(2) S x U(l) (4) 
D SU(2) p x SU(2) S D SU(2) S x U(l) 
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each of which ends in the subgroup SU(2) S x U(l) rep- 
resenting total spin and charge conservation. An SU(4) 
model Hamiltonian containing one and two-body inter- 
actions can then be determined uniquely, 

H = en + vh 2 - G D^D - Givr* • tt 

- xQ Q + kS-S. (5) 

In Eq. (J5J, the parameters Go, G\, and x are effective in- 
teraction strengths of <i-wave singlet pairing, triplet pair- 
ing, and staggered magnetization, respectively, e is the 
average single-electron energy, and v may be interpreted 
as the mean value of the two-body interaction. In later 
discussions, we shall set e = v = because the first two 
terms in Eq. JSJ) are state-independent and provide only 
a constant energy. The term kS ■ S will also be ignored 
for the present discussion, corresponding to assuming the 
total spin for the ground state to be zero. 

Using the properties of Lie algebras, we have shown 
that analytical solutions for matrix elements can be ob- 
tained for the group chains defined in J3J. Each symme- 
try limit has been shown to represent a collective mode 
that corresponds to properties observed in the cuprate 
phase diagram: the SO (4) limit is associated with an AF 
phase occurring when \ is dominant (Go = G\ < x)> 
the SU(2) p limit is associated with a <i-wave SC phase 
occurring when Go is dominant (x = G\ < Go), and the 
SO (5) limit is a critical dynamical symmetry ,103 repre- 
senting a transitional phase that is soft against AF and 
pairing fluctuations (occurring when Go = x). A more 
extensive discussion of these symmetry limits and their 
corresponding phases can be found in Refs. 0, ^| . 

C. Doping dependence in the SU(4) model 

The phenomenology of the cuprate superconductors 
suggests that the expectation value of the system Hamil- 
tonian should depend microscopically on the amount of 
hole doping. Within the SU(4) model it is doping that 
drives the system from one dynamical symmetry limit to 
another, causing the system to undergo crossovers and 
phase transitions. Let us now explore this microscopic 
doping dependence of the SU(4) model in more detail. 

The lowest-order Casimir operators G g for each sub- 
group g in (0} are 
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The SU(4) quadratic Casimir operator is 
C SU (4) =£> f -D + n^ -n+Q-Q + S-S + M(M-A), (6) 



which is an invariant in the SU(4) representation space. 

In the general case we may introduce a seniority-like 
quantum number ni, — wf2, which is the number of parti- 
cles in the system that do not couple to D or tt pairs, with 
u being the number density of the unpaired particles. For 
the u = case, we take as a Hilbert space 

\u = 0) = \n x n y n z n d ) 

= (^(^(Trir^no), 

which is a collective subspace (the D-tt pair space) as- 
sociated with SO(6) irreps of the form (oi, 02, 03) = 
(f ,0,0). [Because SO(6) and SU(4) have the same Lie 
algebra we choose to label SU(4) irreps with SO(6) quan- 
tum numbers.] The expectation value of the SU(4) Cas- 
mir operator JBJ for u = is 

and is a constant for any state in the D—tt pair space. 
More generally, states with rij unpaired particles can be 
described by irreps of the form (|(fi — «&), 0, 0). Thus 
for the w^0 case with broken pairs the expectation value 
of the SU(4) Casimir operator is 

HCso(4) ,„, = §(§ + 4)-|(,-f + 4). 

For a system with charge (M) = |(n — O), the SU(4) 
invariance leads to a conserved quantity 

£ SU 4 = £>*£> + # -tt + Q-Q + S-S 
= G SU 4-M(M-4). 

The expectation value (Ssui) is thus a constant, inde- 
pendent of how the system changes its state within the 
SU(4) space. It follows that in the large fi limit, 




where n is the electron number. In the above expressions, 
x may be regarded as the relative doping fraction in our 
theory. Since il — n is the hole number when n < Q, 
positive x represents the case of hole-doping, with x = 
corresponding to half filling (no doping) and x = 1 to 
maximal hole-doping. Negative x (n > Q) is defined nat- 
urally as the relative doping fraction for electron-doping. 

The true doping rate, defined as P — (Q—n) /fl e , where 
Q, e is the number of lattice sites, is related to x through 



with P > for hole doping and P < for electron doping. 
Pf can be regarded as the maximum possible value of the 
true doping rate P. Experimentally, Pf is found to be 
0.23 - 0.27 0. 
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D. No-double-occupancy and maximum doping 

The SU(4) group is the minimal symmetry accom- 
modating superconductivity and antiferromagnetism in 
cuprates. It is further found that the SU(4) symme- 
try is a consequence of non-double-occupancy - the con- 
straint that each lattice site cannot have more than one 
valence electron. It has been demonstrated [uj that 
with no approximation to the d-wave formfactor g(k), 
the momentum-space operator sets Q and may be 
expressed in the coordinate space as 
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In Eq. JHJ, the summation is over even sites of the lattice, 
the quantity Qij is defined by 



Qij — Qi 



2 u ij 2 & ) 



and j (c r ,i) creates (annihilates) an electron of spin 

i located at r, while c\ i (cf,i) creates (annihilates) an 
electron of spin i at its four neighboring sites, r ± a and 
r ± b, with equal probabilities (a and b are the crystal 
constants along the x and y directions, respectively, on 
the copper-oxide plane), 



Beyond this doping fraction the no-double-occupancy 
condition cannot be ensured and exact SU(4) symme- 
try cannot be preserved. The empirical maximum dop- 
ing fraction 0.23 ~ 0.27 for cuprate superconductivity 
may then be taken as indirect evidence for a strongly- 
realized SU(4) symmetry underlying the superconductiv- 
ity in cuprates. 



III. GAP EQUATIONS IN THE SU(4) MODEL 

Equation defines three dynamical symmetry group 
chains, each representing a collective mode. Analyti- 
cal solutions associated with these symmetry limits have 
been derived However, a realistic physical system 
may not lie in any of these dynamical symmetry lim- 
its. In order to study properties of a realistic system and 
phase transitions between the symmetry limits, one needs 
systematic ways to deal with approximate symmetries. 
There is a well-developed theoretical approach to relat- 
ing a many-body algebraic theory to an approximation of 
that theory that exhibits spontaneous symmetr y br eak- 
ing: the method of generalized coherent states (ToL Il7| . 
This method is a variational procedure using the coherent 
state as the trial wavefunction in a quasiparticle space. 
It may be viewed as the most general Hartree-Fock- 
Bogoliubov variational formulation, but with an addi- 
tional proviso that the variational states are constrained 
to respect the highest symmetry of a set of group chains, 
as in Eq. J3J|. Generalized coherent states are thus par- 
ticularly suitable for studying the ground state proper- 
ties of any strongly correlated many-body system that is 
amenable to a dynamical symmetry description. 
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Explicit commutation operations show that only if the 
no-double-occupancy constraint is imposed, so that the 
anticommutator relation 



{ C M ,4',i'} 



S rr f Si'. 



is valid, does the operator set (JSJ close under the 
U(l) x SU(4) Lie algebra. Thus, the coordinate-space 
commutation algebra of the operators |JS} demonstrates 
that SU(4) symmetry necessarily implies a no-double- 
occupancy constraint in the copper oxide conducting 
plane. This suggests a fundamental relationship between 
SU(4) symmetry and Mott-insulator normal states at half 
filling for cuprate superconductors. 

One immediate conclusion [Tl| is that the no-double- 
occupancy constraint sets an upper limit for the num- 
ber of doped holes if SU(4) is to be preserved exactly: 
!] max = \£l e - Thus, the maximum doping fraction con- 



sistent with SU(4) symmetry is 

ft 
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A. The generalized coherent-state method 

The coherent state \ip) associated with the SU(4) sym- 
metry can be written formally as 



1^)^10*), 
with the operator T defined by 

T = cxp(?7oopI 2 + mWi2 ~ n - c -)- 



(9) 



(10) 



In Eq. |0*) is the physical vacuum (the ground 

state of the system), the real parameters 7700 and 7710 are 
symmetry-constrained variational parameters, and h. c. 
stands for the Hermitian conjugate. Since the variational 
parameters weight the elementary excitation operators 
p\ 2 and q\ 2 in Eq. 110|) . they represent collective state 
parameters for a D-ir pair subspace truncated under the 
SU(4) symmetry (IsJ . 

The symmetry-constrained variational Hamiltonian is 



H' = H — An, 
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where H is the Hamiltonian (JSJ) and A is the chemical 
potential, determined by requiring particle-number con- 
servation. The parameters 7700 and r?io in Eq. ifTHjl are 
determined by the variational principle S(H') = 0, where 
(H 1 ) is the expectation value of H' with respect to the 
ground state |0*) 

(H') = (0*|if'|0*). 

As shown in Appendix A, it is convenient to evaluate 
the variation 5(H') = using a 4-dimensional matrix 
representation that was introduced in Refs. OH Gi|. In 
this representation the unitary operator T implements a 
transformation from the original particle basis to a quasi- 
particle basis and the variation parameters ?yoo and 7710 
are replaced, respectively, by u± and v±, with a unitary 
condition u± + v± = 1. Under this transformation the 
basic fermion operators 

{ c rT> c rJ.) c fh c fl} 
are converted to quasiparticle operators 

{4 T ,a^,a fT ,a f J, 

as shown explicitly in Eq. (| A4|) . This implies that 

(u + C r f + I *) = a r]W 

(u-Cri - W-c| T J |0*) = a rl \ijj) 

(«-<4t + V -°rl) |0*> = 4tI^) 

(w+4 ; - u+c rT j |o*) = 

and one sees that this is a Bogoliubov-type transforma- 
tion: each quasiparticle state is a mixture of a parti- 
cle and a hole, and the coherent state \ip) is an SU(4)- 
symmetry constrained quasiparticle vacuum. 

By using the matrix representation one can calculate 
expectation value for any operator O in the coherent state 
representation through the transformation 

(0*|O|0*) = (^TOT" 1 ^). 

Detailed derivations are given in Appendix A. 

B. Temperature dependence 

At finite temperature, may no longer be a quasi- 
particle vacuum state and the quasiparticle annihilation 
operators acting on \ip) do not necessarily give zero. In 
Appendix B a formalism is derived to deal with the 
finite-temperature case. To formulate the simplest ini- 
tial model we assume that at a temperature T the single- 
particle levels e r ± (defined in Eq. (fTTj) below) are degen- 
erate and contain n r+ + rv_ quasiparticles. The quasi- 
particle number densities are then assumed to be given 



by the Fermi-Dirac distribution 

M t) = £ £ ^ T )= 1 + exp(l ±/fcB rr C n > 

r— even 

where e± is the quasiparticle energy defined in Eq. 
below. In Eq. we have introduced an energy scal- 

ing factor < R < 1. This is because the actual single- 
particle energies are generally non-degenerate so that the 
realistic quasiparticle excitation should be easier. Thus 
n± (T) should in general be larger than those with the de- 
generate approximation. The factor R accounts for this 
effect in an average manner, and it may be determined 
by fitting to data. 

For one-body operators at finite temperatures the ex- 
pectation values are (see Appendix B) 

<£>t) = (D) = -q [P + (T)u + v + + P-(T)u-V-] 

(4) = (ir z ) = -%[P + (T)u + v+-P-(T)u^] 

(Qz) = %[P + (T)v 2 + -P-(T)v 2 _] (12) 

(h) = % [P + (T)(2v 2 + -l) + P_(T)(2v 2 _-l) + 2] 

fa.) = (%> = 0) = (fix) = (Qy) = 0, 
where we have defined 

P ±( T ) = l- M T)=tanh(||). 

For the scalar products of these one-body operators the 
expectation values are products of the corresponding one- 
body ones in the large-Si approximation, 

(D*D) = (D) 2 

<tt + • tt) = (ir z ) 2 (13) 
(Q-Q) = (Qz) 2 
(S-S) = 0. 

If T -> 0, then P± (T) - > 1 and Eqs. O and (JT2J) reduce 
to Eqs. HA9f) and l|A10|) . respectively. 

C. Energy gaps and gap equations 

We now use the preceding results to express the vari- 
ational Hamiltonian (H 1 ) in the coherent state represen- 
tation. Introducing the energy gaps 

A d = G oy J(DW) 

A, = Gi^Trt -t?) (14) 

one obtains 

/ A 2 , A 2 A 2 \ 
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Variation of (H 1 ) with respect to u± or v± (that is, solv- 
ing S(H') = 0) yields 

2u ± v ± {e ± - A) - A ± (n| - 4) = 0, 

which is satisfied by 



in 



e± — A 

e± 

e± — A 
e± 



(15) 



where 



and 



e± = V( £ ±- A ) 2 + A J 



A, 



(16) 



(17) 



Inserting Eq. I|15(l into Eqs. Q12I - I13|l and employing 
the gap definitions 1)140. one obtains the temperature- 
dependent gap equations 



(18a) 
(18b) 



A d = -^-(w + A + + w-A_) 

A, = -^( W+ A+- W _ A_) 
4A 

= W+ (A 9 + A')+^-(A g -A') (18c) 
-2x = w + (A q + X') -w_(A g - A'), (18d) 



where we define 



and 



A' = A - s 



(19) 



(20) 



By solving the above algebraic equations, all the gaps 
and the chemical potential A' can be obtained. The total 
energy can be calculated as 

E = (H') + Xn 



en - 



Ai + A| + Aj 
Go Gi x 



To simplify the discussion, we shall hereafter ignore the 
single-particle energy in the above equation by setting 
e = 0, since this term has been approximated as a state- 
independent constant in our model and thus plays no 
role in the phase competition. We can then express the 
energy density E/Q, as 



E 
Q 



A d 2 , A w 2 A 2 



G n Gift 



=77 ■ ( 21 ) 



The three gaps A^, A^ and A q in the above equations 
represent, respectively, the characteristic energy scales of 
spin-singlet pairing, triplet pairing, and the AF correla- 
tion. Hence the ground state energy is determined by 
the three energy gaps. Once the gaps and the chemical 
potential A' are known, the quasiparticle energies e±, as 
well as the amplitudes u± and v±, can all be determined 
through Eqs. lfT31 - IT7JI . permitting other ground state 
properties to be calculated. 

These results are formally analogous to those of the 
BCS theory with v\ the probability of single particle lev- 
els £± being occupied, A± the energy gaps, and e± the 
quasiparticle energies. The essential difference from BCS 
theory is that conventional pairing theories deal with one 
pairing gap and one kind of quasiparticle; here we have 
two kinds of quasiparticles and several energy gaps, im- 
plying a large variety of new physics. 

In the formalism describing this more sophisticated 
pairing the quantities e± are energies for two kinds of 
quasiparticle excitation, corresponding to two sets of 
non-degenerate single particle energy spectra {e± } sepa- 
rated by an energy 2A g . Each level can be occupied by 
only one electron of either up or down spin. The corre- 
sponding pairing gaps are A±, which are linear combina- 
tions of the two gaps A^ and A w . The probabilities for 
single-particle levels to be occupied or unoccupied are v± 
and respectively. 

In the following sections we shall give analytical solu- 
tions for the gap equations (|18fl , first for zero temperature 
and then for finite temperatures. As we shall see, a rich 
phase structure emerges naturally in these solutions as a 
consequence of competition between the various energy 
scales. 



IV. SOLUTION OF GAP EQUATIONS AND 
THE GAP DIAGRAM AT T = 



There are three parameters, Xi Go, and G±, in the 
coupled algebraic equations l(T%|) . corresponding to the 
three elementary interactions in the SU(4) model: the 
AF correlation (x), the spin-singlet pairing (Go), and 
the spin-triplet pairing (Gi). Physical solutions of the 
gap equations depend on the choices for these param- 
eters. Experimental evidence suggests that these three 
interactions in cuprates are all attractive, and we will 
demonstrate later that the AF correlation should be the 
strongest and the spin-triplet pairing the weakest. Thus, 
in the results presented below we assume 



X > Gq > Gi > 0. 



(22) 



Analytical solutions for the gap equations assuming this 
condition can be obtained as follows. 
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A. Solution of the gap equations at T=0 



(4) The pure AF solution: 



The gap solutions at T — can be written explicitly for 
two doping regimes separated by a special doping point 
given by 



X-Gi 



(23) 



We shall interpret critical doping marking a 

quantum phase transition because the wavefunctions and 
physical properties of the two doping regions lying on 
either size of this point at zero temperature are qual- 
itatively different. Specifically, one finds the following 
solutions. (The general derivations of the gap solutions 
are given in Appendix C.) 

(1) The all-gap solution for x < x q : 



\ = 
A d = 

A. = 
A' = 



G n 



\J x(x q 1 - x) 

\J x(x q - x) 



X q [L X q X) 



Gift 



(24a) 
(24b) 
(24c) 
(24d) 



This is the most important solution that has all gaps 
non-zero and exists only for the doping range x < x q . In 
addition, there are trivial solutions in which at least two 
gaps are zero. 

(2) The pure spin-singlet pairing solution with A q = 
A w = is valid for the entire physical doping range < 
x < 1 and given by 



A 9 
A 



A d = 





Go ft 



A' = - 



Gnft 



(25a) 
(25b) 

(25c) 



This is also an important solution in that it gives the 
ground-state for x > x q (see discussions below). It can 
be verified easily that both solutions (|24|l and (|25|l satisfy 
the SU(4) condition J7J with u = 0, which means that 
all electrons are paired at T = 0. 

For completeness, we list below additional trivial so- 
lutions. These are a spin-triplet pairing solution, a pure 
AF solution, and a metal solution (with all gaps zero), 
all of which are valid for the entire physical doping range 
< x < 1: 

(3) The spin-triplet pairing solution: 



A, = 

Ajf = 

A' = 



A d = 
Gift 
2 

Gift 



(26a) 
(26b) 

(26c) 



A d = A„ = 

xft 



A 



— (!-*)■ 



(5) The metal solution: 



A, = A d = A w = 0, 
A' = -2kT atanh(x) 
= at T = 0. 



(27a) 
(27b) 

(27c) 



(28a) 
(28b) 



The trivial solutions can be verified easily. For exam- 
ple, when A q = one obtains immediately w + — w- 
from Eq. (|18cl) . There are two solutions satisfying Eqs. 
(|T%a|> and (fT8E|) : either A d ^ and A^ = 0, or A T ^ 
and A d = 0. The former gives the singlet pairing so- 
lution l|25al - I25b|l while the latter leads to the triplet 
pairing solution Ij26al - l26bjl . Eqs. Ij25cjl and Ij26cjl for 
A' result from Eq. I|18d|) . The AF solution is obtained 
because A^ = A^ = is a trivial solution of Eqs. Q18al - 
I18b|l . The metallic state solution with all gaps vanishing 
is obvious from inspection of the equations. 

We shall demonstrate below that these solutions con- 
tain rather rich physics as a function of doping. Among 
the five sets of gap solutions (Eqs. I|24l - l2*S|l ). the one 
with the lowest energy at each doping corresponds to 
the physical ground state. We can calculate these ener- 
gies by inserting the gap solutions directly into Eq. il^TJl , 
and then investigate how these different sets of solutions 
compete with each other at T — 0. From Fig.^ on e sees 
clearly that the energy of the all-gap solution E m i x is al- 
ways the lowest one and thus is the physical ground state 
for the doping range x < x q . For x > x q , the pure singlet 
pairing state becomes the ground state because is the 
lowest energy for this doping range. All the other triv- 
ial solutions lie higher in energy. They may be regarded 
as collective excited states but they cannot become the 
physical ground state at T = 0. 

As we will see later in the study of phase transitions, 
for T > the AF or the metallic state could become the 
ground state in certain temperature and doping ranges. 
However, note that this can never happen for the spin- 
triplet pairing state, as long as Gi is the weakest of 
the three coupling parameters. In other words, although 
spin-triplet pairing plays an important role in the SU(4) 
theory, the pure spin-triplet state can never be reached 
by thermal excitations as long as the condition (|22|) is 
satisfied. 



B. Gap diagram at T=0 

A generic gap diagram at T = describing features 
of the energy gaps as functions of doping x is shown in 
Fig. [3 The diagram is constructed using the analytical 
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FIG. 1: (Color online) Total energy associated with different 
gap solutions at T = 0. -E" m i x is the energy calculated with the 
all-gap solution (J24j, while E^ (calculated from the solution 
in Eqs. {25) 1. E^ (Eqs. {5jJ) and E AF (Eqs. {57)1 represent, 
respectively, the energy density of the spin-singlet pairing, the 
spin-triplet pairing, and the AF solutions. The energy of the 
metallic solution is set to be zero and taken as the energy 
reference. The interaction strengths used in this plot are the 
same as those in Fig. |5J 



expressions {24) and {25(1 . It is generic because the rel- 
ative sizes of the gaps can be modified by the choice of 
interaction strengths but their basic forms are dictated 
entirely by the algebraic structure of the model [and that, 
in turn, is determined by the physically-motivated choice 
of generators given in Eq. ©]. Four doping-dependent 
energy gaps are predicted: 

L The gap A q (Eq. {24a) 1 measuring antiferromag- 
netic correlations; A q has its maximal value at 
x = 0, decreases nearly linearly to the region of the 
pairing gaps as doping increases, crosses the pairing 
gaps, and vanishes eventually at the critical doping 



2. The spin-singlet pairing gap (Eq. {24b) 1. which 
is the superconducting gap for x < x q . 



3. The spin-singlet pairing gap Ao = A^ (Eq. {25b) 1. 
which is the superconducting gap for x > x q but is 
not the order parameter for the ground state in the 
doping range x < x q . 
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4. The spin-triplet pairing gap A^ (Eq. {24c) 1. Sim- 
ilar to the case for A q , the triplet gap A„. exists 
only in the doping range x < x q . It has its maxi- 



FIG. 2: (Color online) A generic gap diagram for energy gaps 
versus doping, as predicted by the SU(4) model at T = 0. The 
energy gaps are scaled by Grjfi and the doping parameter is 
scaled by the maximum doping Pf (We assume Pf = 1/4 |Tl|l. 
The interaction strengths are assumed to be x = 13, Go = 8.2, 
and Gi = 1.3 (in an arbitrary energy unit), which, according 
to Eq. 12311 . requires the critical doping point to be x q — 0.64. 



mal value at x q /2 and vanishes at both ends of its 
range. 

Note that the spin-singlet pairing gap exhibits fundamen- 
tally different behavior on the left and right sides of x q : 
for x > x q , it corresponds to a monotonic curve, labelled 
as Ao- However, as the doping decreases from x q the 
spin-singlet gap splits into two curves (labelled as A^ 
and Ao) having very different doping dependence. More- 
over, the splitting correlates strongly with the amount 
of hole doping: the lower the hole doping, the larger the 
splitting. 

The qualitative features of the SU(4) gap diagram 
seem to agree with a large body of recent observations 
in cuprates 0- In particular, the appearance of a crit- 
ical doping point and the splitting of the pairing gap 
in the underdoped compounds are basic predictions of 
the SU(4) model that have some experimental support 



1. Critical doping point and pairing gap splitting 

The occurrence of a critical doping point and the split- 
ting of the SC pairing gap are understood in the SU(4) 
model as a direct consequence of competing SC pairing 
and AF correlation in the doping range below x q . When 



9 



doping is small the AF correlation dominates the SC pair- 
ing. In this case, a state with large AF correlations and 
suppressed pairing is favored in energy, and thus can be- 
come the ground state. Therefore, the SC gap A^ for 
the ground state is smaller and the larger pairing gap Ao 
is associated with an excited state. However, as doping 
increases the pairing correlation grows quickly and even- 
tually dominates. The critical doping point x q is just the 
doping fraction at which the AF correlation is completely 
suppressed. 

The critical doping point x q defines a natural bound- 
ary between overdoped and underdoped regimes having 
qualitatively different wavefunctions. It corresponds to 
the doping point where AF correlations vanish and sep- 
arates a doping regime characterized by weak supercon- 
ductivity and reduced pair condensation energy from a 
doping regime characterized by strong superconductiv- 
ity and maximal pair condensation. The optimal doping 
point (corresponding to the maximum of the SC pair- 
ing gap) has been used extensively in the literature to 
mark the boundary between underdoped and overdoped 
superconductors. It is the doping point where the com- 
petition between the AF and SC correlations leads to 
the maximal T c , but our results indicate that it does not 
mark a boundary between qualitatively different physical 
regimes characterized by qualitatively different wavefunc- 
tions. 

According to the SU(4) model, the experimentally ob- 
served location of the critical doping point can set a 
strict constraint on the relative strengths of the three 
elementary interaction strengths, x? Go, and G±. As 
one can see from Eq. (|23|) . \ must be either greater or 
less than the other two strengths Go and G\\ otherwise 
(x — Gq)/(x — G\) would become negative and no x q 
could occur in the physical range < x < 1. Observa- 
tions require that \ be greater than both Go and G\ in 
order for the normal state at the half-filling to be an AF 
state (Mott insulator). It then follows that Go should 
be intermediate in strength between x and G\ ; otherwise 
(x — Go)/(x — Gi) > 1 will lead to x q > 1, which is 
outside the physical doping range. Thus the SU(4) sym- 
metry and the basic experimental observations require 
that the three interaction strengths must satisfy the con- 
dition x > Gq > G-y. This is the condition that we have 
assumed in (f22J - 

It is instructive to examine two extreme cases of Eq. 
(J2SJ: one is x q = I if Go = G\; the other is x q = 
if x = Go. For these two cases no critical doping point 
x q exists within the physical doping range. These cases 
correspond, respectively, to the SO (4) and SO (5) sym- 
metry limits of the SU(4) model, as has been mentioned 
in Section II. A. For the former case (x q — 1), there is 
no room for the SC phase at any x, suggesting that the 
system in the entire doping range is in the AF phase. 
For the latter case (x q — 0), the system is in the SO (5) 
limit. Although the energy minimum in this limit is at 
A q = 0, the same as for the SC case, it is very shallow 
and the system has large AF and pairing fluctuations. 



The SO (5) energy surface at x = is in fact completely 
flat, develops an energy minimum with A g = as the 
doping x increases, and gradually evolves to an SC-like 
state at x — 1. This behavior is characteristic of a criti- 
cal dynamical symmetry, as discussed extensively in Ref. 



2. Nature of the pseudogap 

The observation of pseudogap behavior in cuprates 
0, has motivated a variety of theoretical discussions 
but there is little agreement on the source of this be- 
havior. Two alternative classes of explanation for the 
origin of pseudogaps have received considerable atten- 
tion 0: (1) The preformed pair picture |22| . which sug- 
gests that the formation of pairs and the condensation 
of those pairs into a state with long-range order happen 
on two different energy scales, with the pseudogap being 
the energy needed to form pairs before pairing conden- 
sation. If the pseudogap originates in preformed pairs, 
it should merge with the pairing gap in the overdoped 
region. (2) The competing- order picture |2jj, which sug- 
gests that the pseudogap is an energy scale associated 
with a form of order that competes with superconduc- 
tivity. If the pseudogap originates in competing order, 
it should not merge with the pairing gap but should in- 
stead drop to zero at a critical doping point where the 
competing order is completely suppressed relative to the 
superconducting order. 

In the solution of the SU(4) gap equations, the result- 
ing energy gaps (scales) are A^ (as well as A ), A ff , and 
A q . Generally the gaps A a and A , may be interpreted 
as SC gaps and the gaps A w and A q may be interpreted 
as pseudogaps. However, our results suggest (as men- 
tioned in Section IV. A and discussed further below) that 
the pure spin-triplet state can neither become a physical 
ground state at T = 0, nor be reached by thermal exci- 
tations. The A w gap can coexist with A q and A^ only 
in underdoped compounds and cannot survive above T c . 
Therefore, within the SU(4) model only the AF gap A q 
is a candidate for the observed pseudogap. 

In the competing order picture, the energy scale that 
competes with superconductivity in cuprates is often 
identified with antiferromagnetism |20j . and this energy 
scale may be called a pseudogap. The pseudogap energy 
scale vanishes at the critical doping, which is indepen- 
dent of the disappearance of the SC gap at P « 0.25. 
Furthermore, experiments have shown (for example, Ref. 
[20|) that the observed pseudogap is intimately connected 
with the AF correlations, which disappear exactly at the 
same critical doping point. The AF gap A q has precisely 
these properties and thus may be interpreted as a pseu- 
dogap. Specifically, the scale A g in the SU(4) model is 
the energy per electron required to break the AF cor- 
relation. Therefore, the nature of the pseudogap in the 
SU(4) interpretation is consistent with the competing or- 
der picture. 
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In the SU(4) model the spin-triplet pairing plays a me- 
diating role in the AF-SC competition. As can be seen 
from Fig. 2, its order parameter A ff is zero at x = 
where the AF correlations dominate, increases as dop- 
ing increases to a maximum at half of the critical dop- 
ing value where the AF and SC correlations are compet- 
ing strongly, and finally disappears at the critical doping 
point x q where the AF-SC competition ends. 

On the other hand, the SU(4) model does not con- 
tradict the general preformed pair picture conceptually. 
If we define preformed pairs as pairs that are formed 
before pairing condensation, pairs with only AF corre- 
lations may be interpreted as preformed pairs. (These 
AF pairs correspond to a mixture of D and 7r pairs in 
our Hilbert space.) This conclusion will become clear 
from the discussions later in Section V. Thus A q may be 
viewed also as the scale of stabilization energy associated 
with preformed pairs that condense into a pure SC pair 
state only after the AF correlations are completely sup- 
pressed by increasing hole doping. A pseudogap arising 
from this preformed pair picture can exist only in the 
underdoped regime since A q vanishes at the critical dop- 
ing point. For overdoped systems (x > x q ) there is no 
pseudogap and thus there are no preformed pairs in the 
gap solutions found here. 



V. SOLUTION OF GAP EQUATIONS AT 
FINITE TEMPERATURES 

The gap equations for T > differ from those at T = 
in that the terms w± (see Eq. acquire a temperature 
dependence 

P±(T) tanh(i? e± /2fc B T) 

w± = = . (29) 

e± e± 

For finite temperature the gap equations could have a 
variety of solutions, even for a fixed doping x. Which 
solution should correspond to the physical ground state 
depends on temperature and doping, and is determined 
by minimizing the energy. The solutions that we have de- 
rived in Appendix C are obtained under the (physically 
motivated) condition x > Go > G\ > 0. Under other 
conditions the gap equations could have different solu- 
tions and one has to solve Eqs. I|18I -I2U ) | for individual 
cases. 

By comparing the energy densities expressed in Eq. 
1)21(1. one can then find the solutions having the lowest 
energy and determine the physical ground-state solution 
for given temperature T and doping x. We examine sev- 
eral cases below. 



3. Scaling property of the gaps 



A. The A d + A„ + A T case 



It is well known that pairing gaps in cuprates scale with 
(T c ) max (the maximum superconducting transition tem- 
perature) for all high-T c compounds studied so far [28l |. 
According to Ref. |20j , there is also strong experimental 
evidence indicating that the doping value of the puta- 
tive quantum critical point is universal in the hole-doped 
cuprates. 

In Fig. |3 all the energy gaps are scaled by GqVL. From 
Eq. 124b|) . we know that for any given Go the gap A^ 
depends only on x q . Therefore, if the critical doping 
point has a universal value the doping dependence of A^ 
should also be universal, since we shall show in the next 
section that (T c ) max is proportional to Go- The gap Ao 
also has this scaling property because, according to Eq. 
(|25b|l . A in the whole doping range depends only on Go- 

However, the other two gaps A q and A v do not scale 
in this way because they depend on x/Gq and G\/Gq, 
respectively. Changing the ratio of x/G and G\/Gq can 
change the size of these gaps; thus the doping dependence 
of A q and A ff for different compounds generally could be 
different. Only if the strength of triplet pairing G\ is zero 
(thus Att = 0), can A q be scaled because in this case 

X/G = 1/(1 - x 2 q ). 

Hence, the scaling property of A q may be taken as an 
indicator of triplet pairing strength in cuprates. Recent 
A q data |20l |24| seem to support such a scaling property, 
at least approximately, implying that G\/Gq should be 
small. 



This case corresponds to the all-gap solution found at 
T = for the doping range < x < x q . For T > 0, 
the temperature and doping dependent gap solutions are 
derived in Appendix C, with the results 



A q = ^rJ {xq 1 -y){x q -y) ~ 
z y 

GoO 
~2 



(30a) 
(30b) 



A d = -^\jx{x q x -y)g{y) 

G i 

A^r = -^-^Jx{x q ~y)g{y) (30c) 

v = -fc^(£- v )-2i2*(30d) 

where y and g{y) are defined through 

x 



y = 



g(y) 



y /I + (T)+I-(T)T(y) 

fx | I + {T)-{x~/yJ* 
V y 2x{x q - y) 



with 



r(y) = 

I±(T) 



y/(x q -y){x q 1 -y) 
Pl{T)±Pf{T) 
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and x q = (x 1 + x q )/2. It can be verified that these Note that in the present case, 



solutions satisfy the SU(4) invariant 



Q 2 



(Ssui) = — [I+(T) - x 2 } 



(31) 



With increasing temperature the system begins to break 
pairs, with the number density of unpaired particles at 
temperature T given by 



1- VMT) = 1 



P 2 (T) + P 2 (T) 



(32) 



More specifically, to obtain the energy gaps for given 
doping x and temperature T, one may adopt the follow- 
ing procedure. From Eqs. (|29|) and i|16|) one obtains 



t _ R^(A q ±\') 2 + A ± < 
2/cb atanh(ui±e±) 
This equation implies that 



(A 9 + A') 2 + A + 



(A,-A') 2 + A_ 



atanh(u> + e + ) 
By solving Eqs. (|18cl 



atanh(w_e_) 
I18dl) directly, one gets 

2A„ 



(33) 



(34) 



w± 



A 9 ±A' 



Now by using Eq. I|C2|) in Appendix C to convert A T into 
A d 



Ajr = 



UJQ 



A d , 



A + and A_ can be related to A a 



A+ = 



1± 



- w 



■ Uli 



(35) 



(36) 



Thus, for a given doping x, one can solve for y from Eq. 



aj, and A' from Eq. (|30djl . for each A q . The pair- 
ing gap Ad can then be obtained by solving Eq. I|34|) 
directly. With A^ determined, A n and the correspond- 
ing temperature T can be calculated from Eqs. (|35|l and 
|J33}, respectively. 



B. The A d case 

This case corresponds to one of the trivial solutions at 
T = 0. Following a similar procedure as in the T = 
case, one obtains the temperature and doping dependent 
gap solutions 



A d 
A' 



A w = 

Go n 

2 

G n 



VU(T)-x 2 



(37a) 
(37b) 

(37c) 



I+(T) = Pl(T) = P 2 {T). 

By using Eqs. dJ, (HJ}, and (|37c|) . Eq. l|37b|) may be 

written in the following form 



J(^x 



T = 



2fce atanh 



( 2A d 
Won 



(38) 



Therefore, for given x and T the gap A^ can be obtained 
from Eq. and A' obtained from (|37c|) . 

Since now all the correlations are zero except the spin- 
singlet pairing, we have 



SU4/ 



(£> f i?) 



Ad 
Go 



One can see that the SU(4) invariant l|31|l is also valid 
in this case, and the number density of unpaired parti- 
cles can be evaluated using Eq. (|32|) . It can be checked 
without difficulty that when T — > 0, 

p±ct)->i j+cn-i, 

and the solutions l|37|) then reduce to Eq. (|25|l of the 
T = solution. 



C. The A 9 case 

In this case both A^ and A v are zero, corresponding 
to a solution that applies only for temperatures T > T c : 







A, = A-(P_(T)-z). 



(39a) 
(39b) 



These results can be derived from Eqs. H18c|l and IjlSdJl . 
which in the present case reduce to 



4A, 

2x 



P_(T) + P+(T), 
P-(T)-P+(T). 



(40) 
(41) 



Using Eq. (|41|) . Eq. (|39b|l can be rewritten as 



In a manner similar to that for the Ad case, one can show 
that Eqs. (|31l) and (|32|l are also valid for the pure AF case 
(A d = A w = 0). Thus Eqs. (ED an d ((321 are actually 
general expressions; for different cases only the values of 
the quasiparticle energies e± (contained in P±(T)) differ. 
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By solving Eqs. I|40|) and l|41|l one obtains 

RA a 



with 



-4, 



atanh 



T 



A' 



2A g 



k B A + 
k B TA_ 
R 



± atanh 



(42) 
(43) 



/2A„ 



Thus, for a temperature T, A q can be obtained from Eq. 
(gH and A' from Eq. igSJ. 

There is an intriguing point concerning the pure AF 
state at T = that should be clarified. On one hand, 
when T = 0, we have that P_(0) = 1 and P+(0) = 1-22. 
Thus 

7+(0) = l-2iE + 22- 2 , 

and according to Eq. 132|) . u > 0. This means that, if 
i/0, the pure AF phase requires pair breaking even at 
T = 0. The SU(4) invariant in this case is 



SU4/ 



(O f O> 



" 2 n ^ 2 



On the other hand, it is our basic assumption that when 
T = all particles are paired (it = 0) and the SU(4) 
invariant should be 



SU4/ 



4 



which is state-independent if SU(4) symmetry holds. De- 
spite appearances, these two expressions for (Ssui) & re 
not in contradiction because when there exist pairing cor- 
relations (one or both of Go and G\ nonzero) , no AF state 
can be the T — ground state unless x = 0. The only 
possibility for the pure AF state to become the ground 
state at T = for x > is in the absence of pairing 
(Go = G\ =0). But under this condition it is not neces- 
sary to require (D^ D) = (tt^tt) = 0, so that (D^D) and 
(tt^tt) can supplement (Q^Q) to preserve the invariancc 



SU4, 



= f (1 



VI. THE PHASE DIAGRAM 

We now use the preceding results to construct temper- 
ature versus doping phase diagrams in the SU(4) model. 
The interaction strength parameters of the theory rep- 
resent effective interactions within a truncated space. 
For an effective interaction approach to make sense, the 
interaction strengths cannot have very strong local de- 
pendence on global parameters such as doping fraction 
(except for possible rapid changes at phase boundaries). 
However, since the effective interaction strengths repre- 
sent renormalized interactions within a highly truncated 
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FIG. 3: (Color online) Phase diagram predicted by the SU(4) 
model with R — 0.6. The interaction strengths are the same 
as those used in Fig. [5] but in units of fes(T c )max, where 
(Tc)max is taken to be 90 K. The two points marked as a 
and b are tricritical points. The critical doping point is at 
P q = 0.16 (corresponding to x q — 0.64). 



space, they may reasonably be expected to have a smooth 
dependence on the global parameters. Nevertheless, it 
is useful to make the simplest assumption for a start- 
ing point: the effective interaction strengths are constant 
across the entire doping range relevant to cuprates. Al- 
though this is a more stringent restriction than is war- 
ranted for a realistic theory, it has the advantage of sim- 
plicity and it should give at least sensible qualitative re- 
sults if the present approach is valid. 

In Fig. we show a typical phase diagram. In the 
calculations the interaction strengths are kept constant, 
with the same values as those in the zero-temperature 
case discussed in Figs. 1 and 2. The only adjustable pa- 
rameter is R (appearing in Eq. lfTT|l h the energy scaling 
parameter that approximately corrects for the possible 
error caused by the degeneracy assumption in the quasi- 
particle excitation spectrum for a finite T. 

There are four distinct phases emerging in Fig. 
a pure antiferromagnetic phase (AF), a superconduct- 
ing phase (SC), a transitional phase with all three cor- 
relations present, which may be called a mixed phase 
(marked as AF+SC), and a metallic phase. The cor- 
relations (energy gaps) associated with each phase are 
indicated in parentheses. The doping-dependent transi- 
tion temperatures T c , T* , and T q define the boundaries 
for these phases. There are two tricritical points a and 
6, which are the intersection points of three phases. 
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A. Phases below the critical temperature 

We first discuss the phases below the critical temper- 
ature T c . In the phase diagram of Fig. 3 there are two 
phases below T c , the mixed phase and the SC phase. The 
essential difference between them is that the mixed phase 
has both AF correlation A q and triplet pairing A w mixed 
with the singlet pairing A^, while the SC phase has a 
pure Ad correlation with A q = A n — 0. At zero tem- 
perature, the separation point between these two phases 
corresponds to the critical doping point x q . 

The boundary between the mixed phase and the SC 
phase is marked in Fig. 3 as T q , which is the locus of 
all points where the corresponding energies of the two 
phases are equal 

[A^CT,)] 2 = [A^(T g )] 2 + [At 1 ^)] 2 + [AW(T,)]^ 
Go Go Gi x 

(44) 

Here we use the superscript indices '(1)' and '(2)' to dis- 
tinguish quantities evaluated in the mixed and SC phases, 
respectively. 

Eq. I|44l) indicates that the order parameters A^, A,r, 
and A q could be discontinuous across the boundary ex- 
cept for T = 0. This implies that the mixed-SC phase 
transition is second-order at T = 0, but could generally 
be first-order for finite T. The actual situation depends 
on the relative location of the critical points a, b, and x q . 
Using P a (or x a ) with a — a or b to represent doping 
fractions corresponding to the tricritical points, we can 
show that the transition is indeed first-order if x q < Xb, 
as illustrated in Fig. 3. 

Quite interestingly, if x q > Xb the point a merges with 
the point b, as illustrated in Fig. 0] The point b is now 
a quadcritical point, where at this unique doping and 
temperature all four phases that we have identified in the 
SU(4) model can coexist in equilibrium with each other. 
In this case, all the order parameters change smoothly 
across the boundary and the mixed-SC phase transition 
becomes second-order. The curve T q in this case looks 
like an extension of the curve for transition temperature 
T*. The results in Fig. 4 are obtained for a larger x q 
value than for those in Fig. 3. This suggests that the 
precise value of the critical doping point is related to the 
number of tricritical points in the phase diagram and to 
the order of the phase transitions. 

Both results presented in Figs. 3 and 4 are possible so- 
lutions. The only difference between them is the triplet 
pairing strength G\. In Fig. 3, G\ = 1.3 is used, resulting 
in a smaller x q — 0.64, while in Fig. 4, G\ = 4.7 corre- 
sponding to x q — 0.76. The former case has x q < Xb and 
the latter x q > Xb- This suggests that the relative posi- 
tion of x q and Xb could be an indicator of the strength 
of triplet pairing. It would be extremely interesting to 
confirm experimentally whether the real phase diagram 
in cuprates is of the Fig. 3 or Fig. 4 type. 
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FIG. 4: (Color online) Phase diagram predicted by the SU(4) 
model for P q > P^. P q — 0.19 (corresponding to x q — 0.76) 
is chosen in this figure while Pb = 0.175 (corresponding to 
x q — 0.7) is the same as that in Fig. [3] All the parameters 
remain unchanged except that Gi increases from 1.3 to 4.7 
(in units of fcs(Tc)max) for a larger P q value. 



B. The critical temperature 

Let us use Fig.|3|to discuss the critical temperature T c . 
The curve for T c in the SU(4) phase diagram consists of 
three segments: 

( T mix (P<Pa) 
T c = { T dl (P a <P< Pb) , 

{ T d2 (P > P b ) 

which are seen to behave quite differently The discus- 
sions below apply also to Fig.0|because in this case there 
is only one phase boundary intersection point b and thus 
the segment Tdi does not exist. 

In the doping range P > Pb the critical temperature is 
T c = Tc\2, which is the SC-metal transition temperature 
and can be expressed analytically using Eq. (|38|l with 
A d = 0: 

T c = T d2 = T(SC <-> metal) 
Rx 

4/cb atanh (x) 

In the doping range P < P a the critical temperature 
is T c = T mix , which is the mixed-AF transition tem- 
perature determined from Eq. i.YMi with A^ = (thus 
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0): 



T„ 



T(mixed <-> AF) 



i?|A,±A'| 



/2A 
2fce atanh ( — ^ =p a; 



(a; < x a ), (46) 



where A' depends on A q through y [see Eqs. (|30a() and 
(|30djl ] and A q is determined by Eq. ipHI) . 

In the doping range P a < P < Pb the critical temper- 
ature is T c = Tdi, which is the SC-AF phase boundary 
and determined by the condition that at T c the energies 
of the two phases are equal: 



[A d (T c )] 2 {A q (T c )f 



Go 



X 



(47) 



Using this condition and Eqs. I|38[> and 142(1. one can 

obtain the pairing gap A d through the following equation 



2 atanh 



'GM 



A, 2 



Vx/GoA d 



, / 2A d \ / 2A d 

atanh | — , — x I + atanh [ — , + x 



and T c is then determined by the equation 
T c = T d i = T(SC «-» AF) 



2 fee atanh 



Y 2A d 



(48) 



which is valid in the doping range x a < x < x\>. 

The phase transitions in the two regions (P < P a 
and P > Pb) are second-order (all gaps change smoothly 
across the boundary). However, transitions in the range 
x a < x < Xb behave very differently. The pairing gap 
Ad in this doping range drops discontinuously to zero 
when crossing the boundary, while the AF gap A q jumps 
from zero to a finite value, suggesting that the system 
undergoes a first-order phase transition from the SC to 
AF phase. 

In traditional BCS theory the energy gap appearing in 
the density of states is identified with the pairing gap. 
This remains true for the overdoped regime (T^) where 
the singlet pairing is the only correlation. However, this 
is not true for underdoped systems where more correla- 
tions are involved and a pairing gap is not necessarily 
equivalent to the energy gap appearing in the density of 



states. As shown in Eq. (|16() . in the underdoped regime 
with the mixed phase there are two gaps, A±, in the 
density of states, but neither of them is an order param- 
eter. The actual process of the mixed-AF transition is 
as follows: with increasing temperature, A n — * at a 
temperature that causes u>_ = wq before T c is reached 
[see Eqs. (|35l - l36fl ]. and at this temperature A± reduce 
to Ad- After that, the order parameter A d — > at T c . 

This complexity in gap structure suggests that one has 
to be careful with the interpretation of energy gaps mea- 
sured in underdoped compounds. According to the SU(4) 
model, there are several gaps involved in this doping and 
temperature range such as A d , Ao, A„, A + , and A_, 
and they all can correspond to energy scales having the 
same order of magnitude. Therefore, an energy gap by 
two different experimental probes may not really be the 
same gap since a particular measurement may be sensi- 
tive only to certain of these gaps because of their micro- 
scopic structure. 

The location of the tricritical point x a can be deter- 
mined from T m i x (a; a ) = Tn(x a ) with the condition that 
the energies of the mixed and SC phases at this point 
are equivalent. By using Eqs. ffify and |@SJ, x a can be 
obtained by solving the following equations 



'G n 



\A q ±X'\ 



+ A d 2 



atanh [ — 



atanh 



( 2A d 

\G Q n 



where A q and A' are known if T mix 



is determined, while 



Ad is calculated through the condition A q /\ — A 2 /Gq. 



C. Phases above the critical temperature 

For temperatures above T c there are two possible 
phases: the pure AF phase and the metallic phase. 

As already mentioned, the pure AF phase at T — can 
exist only at x = as a limit of the mixed phase. This is 
because of pairing correlations that give the mixed phase 
a more favorable energy, and is qualitatively consistent 
with the observation that the AF phase is favored only 
in a very restricted doping region near half filling. 

When T > T c , however, the pure AF phase can become 
the ground state in the underdoped regime, as predicted 
in Fig. |3| This is because thermal fluctuations can de- 
stroy the pairing correlations preferentially relative to the 
AF correlations. If temperature increases further, the AF 
correlation will eventually be totally destroyed (A q — * 0) 
and the system will undergo an AF-metal transition. 

The AF-metal transition is second-order. The tran- 
sition temperature T* is determined by Eq. <|42|l in the 
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limit that A q — > 0: 



T* = T(AF <-> metal) 
RgfLjl - x 2 ) 



(0 < x < x b ). 



(49) 



The Neel temperature can be determined from the pre- 
ceding equation with x = 0, 



Tn = 



4fc R 



(x = 0). 



The location of tricritical point Xb can be determined 
through the condition T^ixb) = T*(xb)- Combining Eq. 
gHJl with gives 



£6=4/1- 



X atanh^t) 



D. Discrepancies 

The phase diagrams of the SU(4) model are consistent 
with much of the current understanding for cuprate sys- 
tems. However, there are two obvious quantitative dis- 
crepancies between the simplified results presented here 
and data. 

One discrepancy occurs at very low doping. Data show 
that the pure AF phase is not confined to half filling 
but can extend over a narrow non-zero doping range 
at low temperature. Experimentally, T c goes to zero 
at P w 0.05, which differs from the theoretical predic- 
tion P = 0. However, recall that we have deliberately 
employed an oversimplified model here (effective inter- 
actions independent of doping) in order to emphasize 
that the qualitative features of the SU(4) quasiparticle 
solution follow from the physics encoded in the algebraic 
structure, not from detailed parameter adjustment. As 
we have shown, in the SU(4) model the condition for the 
pure AF state to be the ground state is the complete 
absence of pairing interactions. Since we have assumed 
all effective interaction strengths to be constant over the 
whole doping range, the pairing correlation exists at all 
dopings. This is not a favorable condition for the AF 
phase to be the ground state and with these assumptions 
T c can go to zero only at the both ends of the doping 
range, x = (P = 0) and x = 1 (P = 0.25). This dis- 
crepancy suggests that there exists an onset of effective 
pairing correlation around Pi ~ 0.05. If we allow a simple 
variation of pairing strength with doping implying that 
the pairs become stable and thus have pairing correla- 
tion only when P > P,;, this discrepancy can be resolved 
easily. 

A second quantitative discrepancy is that our predicted 
Neel temperature is too low (175 K) for a cuprate system 
with (T c ) max sa 90 K. Again, if we relax the assumption 
of constant interaction strength by considering different 
AF correlation strength x before and after the onset of 



pairing near P — 0.05 (for example, maximal x at i = 
that decreases as doping increases and finally is stabi- 
lized at a smaller value when pairing is established), this 
problem is also easy to resolve. Work to establish the 
mechanism for pair formation and thereby to quantify 
the expected doping dependence of the effective interac- 
tions is in progress and details will be reported elsewhere. 



VII. SUMMARY 

The present work has used dynamical symmetries, Lie 
algebras, and generalized coherent states to derive the 
temperature-dependent gap equations expected for a the- 
ory in which antiferromagnetism and d-wave supercon- 
ductivity compete on an equal footing, and for which 
the normal undoped states have Mott insulator charac- 
ter. Although the SU(4) model is constructed using sym- 
metry principles, the coherent state method permits the 
problem to be cast in the form of a generalized quasiparti- 
cle problem. Therefore, the results are presented in terms 
of equations that may be recognized as a generalization 
of the BCS formalism to include more than one kind of 
pairing, subject to an SU(4) symmetry constraint. This 
symmetry constraint has a clear physical origin. As the 
preceding discussion has shown, SU(4) symmetry may be 
viewed as concise shorthand for competing antiferromag- 
netism and d-wave superconductivity on a 2-dimensional 
lattice with no double occupancy. 

The quasiparticle structure that results for the cuprate 
superconductors is rather rich, giving rise to novel gap 
structures including both pairing gaps and pseudogaps, 
and a potentially complex phase diagram. A critical 
doping point P q appears naturally in the theory as the 
boundary between doping regimes having qualitatively 
different ground-state wavefunctions at zero tempera- 
ture: a pure superconducting solution at higher dop- 
ing (P > P q ) and a solution with superconductivity 
strongly suppressed by antiferromagnetism at lower dop- 
ing (P < P q ). Thus, the critical doping point is associ- 
ated with a quantum phase transition. The pairing gap 
has been shown to have two solutions for P < P q : a small 
gap, associated with competition between superconduc- 
tivity and antiferromagnetism that is responsible for the 
ground-state superconductivity in underdoped systems, 
and a large gap without antiferromagnetic suppression 
that corresponds to a collective excited pairing state. 

Within the SU(4) model a pseudogap can occur natu- 
rally. It has been demonstrated that the pseudogap orig- 
inates from the competition of antiferromagnetism with 
e?-wave superconductivity in the underdoped regime, and 
terminates exactly at the critical doping point. These 
conclusions are in accord with many current observations 
in cuprates. Although the pseudogap arises directly from 
competing antiferromagnetic and superconducting order, 
we have also argued that it may be interpreted in a pre- 
formed pair picture since the corresponding wavefunction 
contains singlet and triplet pairs fluctuating in an antifer- 
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romagnetic background but no long-range pairing order. 

Once the parameters in the SU(4) gap diagram are de- 
termined by fitting to the gap data, a rich phase diagram 
as a function of temperature and doping can be sketched 
with only one additional adjustable parameter. A vari- 
ety of phases reflecting the interplay of the spin-singlet 
and spin-triplet pairing with the antiferromagnetic cor- 
relations has been predicted. Properties of the phases 
may be expressed quantitatively using the mathematical 
properties of the SU(4) algebra and its coherent states, 
leading to a set of testable predictions concerning phase 
structure and phase transitions in the cuprate supercon- 
ductors. 

Finally, we note that the striking properties in the 
phase diagram discussed in this paper have similari- 
ties with properties observed in other materials (often 
at lower temperature and energy scales). For exam- 
ple, in heavy-fermion compounds and some 2D organic 
superconductors a superconducting phase appears near 
the boundary of an AF phase. As a second example, 
the manganites have complex competing phases, some 
bearing a resemblance to those discussed in this paper. 
Therefore, we expect that the general formalism pre- 
sented here, which is uniquely suited to deal quantita- 
tively with multiple competing low-temperatures phases 
in strongly correlated systems, will prove applicable to a 
much broader range of systems with doping replaced or 
supplemented by additional control parameters such as 
pressure or strength of magnetic field. 
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APPENDIX A: MATRIX REPRESENTATION OF 
THE SU(4) ALGEBRA 



where O kl (k,l = 1 or 2) are 2x2 matrices with matrix 
elements ojj 5 . A 4-dimensional faithful matrix repre- 
sentation of SU(4) generators © can be obtained imme- 
diately from this mapping: 
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where a x , <r y , and a z are Pauli matrices in the standard 
representation, a± = \{o x H&y), and / is a unit matrix. 
A corresponding matrix representation of the operators 
© is then constructed readily from Eq. I|A2I) . The uni- 
tary transformation operator T of Eq. (jSJ may be written 
in this matrix representation as 



T 



with the definitions 



Yi X 

-Xt Y 2 



(A3) 



v+ ' 




' u+ 




"u_ 


-v- 


Y x = 


u- 


Y 2 = 


u+ 



where the requirement of unitarity implies that 



To facilitate the evaluation of S(H') = in section 
III. A, it is convenient to use a matrix representation 
introduced in Ref. jToL Il9| . In this representation the 
SU(4) generators are expressed in terms of 4x4 matri- 
ces, with the matrix elements defined in the following 
4-dimcnsional single-particle basis: 

{4|0*>, 4,10*), c fT |0*), c n |0*>}. 



Explicitly, for an operator O this is implemented 
through the following mapping: 



Qll Q12 

Q2X Q22 



o 
o 



E 



[o (11) c t. 



(12) + J 
ij r,' rj 



O i21) c- c ■ 



(Al) 



u z ± +v z ± = 1. 

The u's and v's are variational parameters in the matrix 
representation and are related to the parameters ryoo and 
r/io in Eq. (fTOf) . 

The Bogoliubov-type transformation l|A3(l applied to 
the D-tt pair space may be viewed as a quasiparticle 
transformation that is further constrained to preserve 
the SU(4) symmetry. The physical vacuum state |0*) 
is transformed to a quasiparticle vacuum state and 
the basic fermion operators, 



rT> C ri' Cr T' c rl}i 



are converted to quasifermion operators, 
K^a^aft-"?!} witn a ^W 



0, 
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through the transformation 



and for two-body terms 
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/ c rT 
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(A4) 



Using the transformation (|A4() and the mapping QAlj l , 
one can express any one-body operator in the quasipar- 
ticle space as 

0(11) 0(12) 
0(21) 0(22) 



(22) f 

ji a fi a rj 



TOT- 1 = 

r,i 

+ Olfalatj+Oifa*^}, (A5) 

where we put a prime on the summation symbols to indi- 
cate that the summation runs only over r £ even lattice 
sites. The are fixed by the transformation proper- 

ties of the operator O : 



O 



{ixv) 



m.n 



(A6) 



and T^ m ) and are two-dimensional submatrices 

of T and O, respectively. 

Because the quasiparticle annihilation operator acting 
on the quasiparticle vacuum gives zero, the expecta- 
tion values for one-body operators O are given by 

(O) = MdW) = E °u 2) = E Tr(G (22) ), (A7) 

r,i r 

and for two-body operators OaOb-, 

(d A o B ) = {^\6 A 6 B \i>) 



ETr(0f } )ETr(0g 2) ) 



(A8) 
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(21)^(12) 
A B 



Utilizing Eqs. (|A3f) and (|A6|) JA8J), and noting that the 
summation Y.' r provides a factor of 0/2 because the ma- 
trix elements of Eq. (|A6|) do not depend on r, one obtains 
the expectation values for the SU(4) generators and their 
scalar products in the coherent state representation. For 
one-body terms, 
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where we have applied a large-fi approximation, ignor- 
ing terms that are of order l/f2 smaller than the leading 
terms. Under this approximation, the expectation value 
for a two-body operator is simply a product of the ex- 
pectation values of two one-body operators (that is, the 
last term on the right side of Eq. i|A8(l is negligible) . 



APPENDIX B: DERIVATION OF THE 
TEMPERATURE-DEPENDENT GAP 
EQUATIONS 

The formulas in Appendix A are valid only at zero 
temperature. In this Appendix we extend the derivation 
to the finite temperature case. Let us begin with a more 
general Hamiltonian, assuming that the single-particle 
energies are degenerate only for adjacent sites: 



H = E E rhr ~ (GqD^D + GiTT 1, • 7T 
r— even 

+ x Q-Q + nS- 5), 



(Bl) 



where 



n r = 4 T Crf + cljCfl + cl l C Tl + 4 T C fT . 



This Hamiltonian is not an invariant with respect to the 
SU(4) symmetry generated by the operators of Eq. JSJl. 
However, with a slight revision, the following 16 oper- 
ators form an U(4) algebra denoted as U4(r) (with an 
SU4(r) subgroup), 
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As for Eq. ©, one can also define more physical opera- 
tors D(r), 7r(r), Q(r), S(r), M(r), and n r = 2M(r) + 2, 
and express them as 

£> = ]T D(r) t?= J] 

r— even r=even 

Q = ]T S= ^ S(r) 

r— even 



r— even 



The Hamiltonian 1B1I) has the symmetry of the direct 
product of SLUtr). The coherent state now becomes 

\lp) = e (')oo(i-)-D t (r)+r, 1 o(r)7rt(r)-h. c.) | q*^ 

r— even 

= T\0*), 
where T is a direct product of the T(r), 



n nr)= n 
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The variational Hamiltonian is 
H' = H — An 



r— even 

(GqD^D + GxTft • 7? + xQ ■ Q + kS ■ S 



Taking expectation values of operators as in Eqs. (|A9|) 
and HAlOj l permits the expectation value of H' to be 
determined. However, to consider the temperature de- 
pendence, one should replace the quasiparticle vacuum 
state |?/>) by the state \tp(T)) in which quasiparticles are 
thermally excited. The expectation values of a^a r j and 
a li a fj i n <|A5|) are no longer zero when i = j. Therefore, 
instead of (|A7|I , the formula to evaluate expectation val- 
ues for a one-body operator becomes 



(6) = (^(T)|0^(T)) 
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r— even 

[0^\r)h^{T)-0^\r)n n {T)] 

Pg 1] (r)n rl (T) - 0<f (r)rl r - T (T)] } . (B2) 



In Eq. (|B2(1 . we have supposed that at temperature T, 
n r i (T) and hf% (T) are respectively the numbers of quasi- 
particles at r and f with spin i, 

(V(T)|4,a r #(T)) = n„(T) 
(^(T)|o^aH|^(T)) = n fi (T). 

It can be shown that for all the SU(4) generators except 
spin S, 
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where 
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As for Eq. (|A9|I . one can use Eq. (|B3|) to obtain the one- 
body terms 
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where 



P r ±(T) = l-nr±(T). 



As in Eqs. JHI, {n x ), (tt v ), (Q x ), (Q y ), and (S) are 
all zero. Note that for spin this is true only when the 
quasiparticles result from thermal excitations. This is 
because from (IB2I) one can show that 



(S(r)) = (S z (r)) 

which is generally nonzero. Only when the quasiparticles 
are due to thermal excitations, as we will show later, 
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we do have n r -\ = hpi and h r [ — hf-\, which leads to 
(S z (r)) = 0. 

Applying the large-fi approximation (ignoring the sec- 
ond term in Eq. (|A8|) when Q — > oo), one obtains the 
two-body terms 
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= \ E [Pr+(T)u r+ V r+ 
^ r— even 



+P r - (T)Ur-Vr-] 

(#.7?) = (tt z ) 2 



(B6a) 



QQ 



= \ E [P r+ (T)u r+ v rA 

^ r— even 

1 2 

—P r -(T)u r _v r -} > 

= (Q 2 > 2 

= i{ E C p ^( T )(< 

^ r— even 



(B6b) 



M ) = 



-P r _(T)( V 2 _- M 2 _) 



a;fj 

2 

i £ [p, + (t)K 2 + - 



(B6c) 



r— even 

+P r _(T)(v 2 r _ - M 2 _)] 



(B6d) 



The next operations are similar to those of section III.C. 
By introducing the energy gaps 

A d = \G (D)\ A n = \G (7r z )\ A q = \ x (Q z )\, 

defining e r ± and A±, and utilizing (|B5|) to construct (H 1 ) 
and perform the variation calculation 8{H') = 0, one 
obtains 



2u r ±v r ±(e r ± - A) - A ± (u 2 , ± - v 2 r± ) = 0. 
Solving this equation gives 
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with quasiparticle energies defined as 



e r ± = \J (e r ± - A) 2 + A± 5 
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Inserting Eqs. (|B8fl into Eqs. (|B5fl and (|B6|) . one finds 
the gap equations: 
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The energy of the system is then obtained as 
E(T) = (H') + ?iA 
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Utilizing Eqs. (|B8|)- 
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r— even 



This relation can be used to rewrite the energy as a sum 
of the quasiparticle vacuum energy and the quasiparticle 
excitation energy 

E(T) = E {T)+ {n r+ (T) e r+ + n r _(T) e r _ }. 

r— even 

(Bll) 

In the above equation, the quasiparticle vacuum energy 
is 

E Q (T) = Y ( 2s r ~ e r+ -e r _) - (fi- n)A 

r— even 

A d 2 , A, 2 , V 

X 



Go Gi 



(B12) 



T=0 



which is determined by solving the gap equations (|B10|) 
at T — 0, while the quasiparticle number operators are 
assumed to be 

h r ±(T) = 



1 + exp^ 



\k B T) 
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since they are associated with thermal fluctuations, and 
thus 



( er± 
\2k B T 



P r ±(T) = 1 - n r ±{T) = tanh 



Note that, according to the definition 1|B4J| . 

h r+ (T) = n rT (T) + fifi(T) 
h r ^(T) = h r i(T) + hfi(T). 

Therefore, Eq. 1B11|I implies that the quasiparticle states 
\r 1) and |f J.) are degenerate with excitation energy e r +, 
while the states \r [) and |f |) are degenerate with exci- 
tation energy e r _. Thus, as we have mentioned before, 



n r] (T) = n fl (T) 



n rl (T) = n FT (T) 



1 



1 + exp 



e r + 



1 



1 + exp 



e r - 
k B T 



The physical meaning of A,j can now be better under- 
stood. From the expression for the quasiparticle energies 
e r ±, one sees that the original four- fold degeneracy of 
the single-particle energy level e r has been split into a 
two-fold degeneracy because of the Q ■ Q interaction: 



-r+ &r Aq 



(H) and |f|) 
(|ri) and \f\). 



The energy difference 2A g is just the energy required to 
flip the spin of an electron from \r f) to |r J.) or from |f |) 
to \ f |), and vice versa. 

The above results can also be obtained through a mean 
field approximation. Note that any two-body operator of 
the form OaOb can always be written as 

OaOb = [{6 A ) + {6 A -{6 A ))] 

x [(Ob) + (Ob - (6 B ))] 

= (dA)(d B } + o A (dB-(o B )) 
+ (6 A - (6a))Ob 
+ {6a-{6a)){Q b -{6 b )). 

Suppose that the last term of the preceding equation can 
be ignored, implying that the physical effects of (O a — 
(O a}) and (Ob ~ (Ob)) are small. Then the two-body 
operator reduces to an one-body operator: 

A d B = -(0 A )(d B ) +Oa(Ob) + (d A )d B . 
Applying this equation to the Hamiltonian IIBlll . 



H = H + J2 h ( 7 



where Hq is a c-number 



r— oven 



A/ 
Go 



Att 

Gi 



A 



(B13) 



the subscript denotes a ground state expectation value, 
we have made the replacements 

-G (D)^A d , —Gi(ir z ) A,,., x (Qz) A„ 

and the second term in H defines a mean field 

h(r) = 2(e r -X)M(r) + {A d [D\r)+D(r)} 

+A^[7rt(r)+7T z (r)]-2A q Q z } . (B14) 

Now we implement a quasiparticle transformation on H, 



H = THT- 1 = H 



£ Mr), 



where 



h(r) = Th(r)T- 1 = 



U 21 \r) U 22 \r) 



and 



e r+ 
e r - 



h^\r) = 
U 12 Hr) = h^( r y 



-e r _ 
— e r . 



o A 



-O. 







with 



o+ = 



e r ± = 



A±(u 2 r± -v 2 r± ) 
- (e r ± - A) 2u r ±v r ±, 
(e r± -\)(u 2 r± ~v 2 r ±) 
+ A±2u r ±v r ±. 



(B15) 



(B16) 



The condition for the Hamiltonian H to be diagonal is 
that Eq. I|B15() be zero, which is equivalent to Eq. (|B7|1 . 
The solutions of u r ± and v r ± are the same as those in 
Eqs. (|B8p . Inserting these u's and v's into Eq. (|B 16|> . one 
can check that the e r ± appearing in l|B10(l are indeed the 
quasiparticle energies defined in (|B9(1 . 

In the expression for an operator [see Eq. I jAlj l]. h(r) 
can be written as 



h(r) = er+al^cir^ + e r _aj^a r j_ 

— —(e r+ + e r _) + (n r+ e r+ + h r -e r J). 



Since 



H - £(e r+ + e r _) 
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[see <|B12ll and <|B13ll ]. the Hamiltonian may be written 
as 



7t — Eq 



■—even 



77 ^~"T-\- 



(B17) 



which is Eq. (|ljll|) . If we ignore single-particle energy 
splitting, it and v have no r dependence, meaning that 
Y^ r = 0/2, and in this case all results reduce to those 
discussed in Section III. 

From the preceding discussion it becomes clear that 
the SU(4) coherent state is in fact a mean field solution 
of the SU(4) Hamiltonian. The SU(4) coherent state may 
thus be viewed as a generalization of the BCS theory, and 
the process to derive it can be characterized as a Hartree- 
Fock-Bogoliubov (HFB) transformation that implements 
a (non-abelian) symmetry constraint on the variational 
wavefunction. It is a generalized BCS theory because 
the system contains two kinds of quasiparticles instead 
of one. In addition, the Q ■ Q interaction splits the single- 
particle energy levels into two sets, 

e r ± = e r =F A,. 

These new features are primarily responsible for the dif- 
ferences between high-T c superconductivity and conven- 
tional superconductivity described by the ordinary BCS 
theory. 



APPENDIX C: GENERAL SOLUTIONS OF THE 
GAP EQUATIONS 

The temperature-dependent gap equations are coupled 
algebraic equations. In this section we present one way 
to solve these equations. First, we rewrite Eqs. I|18a|) and 
(|18b|l in the form 

(ujq - w+) A d + (oji - w + ) = 
(uj a - W-) A d - (uji - 77J_) A n = 

with 



LOq 



G n 



UJl = 



G x fi 



The condition for these coupled equations to have solu- 
tions is 

(uq — w + ) (oji — 7«_) + (oji — w+) (uj — t/j_) = 0. 
One then obtains 

~0 _ UJq + UJl _ UJqLUi 



W± = U) ■ 



w T — UJ 
— UJ 



UJ 



(CI) 



and 



A 7r =[^^)A d 

W- — UJl 

A-,=-[^±^)A l . 

1/7+ — UJl 



(C2) 



By using Eqs. i|18|) and the new equations IjClf) and l|C2|l . 

a formal solution in terms of w± can be expressed as 



A d = G n 



A, 



(I- 


x 2 )/A 


- (A g /x^) 2 


1 - 




w + /uj \ 




w + /uji J 


(I- 


x 2 )/4 


- (A g /x^) 2 



\ 1 



a = ±n 



1 — W+ 1 UJl 

1 - w + /ui Q 

x{w + — W-) 



A' 



X^lw+W- — (w+ + 7«_) 

xQ(w+ + W-) — 4 



X^W+W- — (w + + 7«_) 



(C3a) 

(C3b) 

(C3c) 
(C3d) 



where we have made use of the SU(4) invariant [see Eq. 
0] 



a 2 



SU4/ 



(I-x 2 ) I={l-uf 



(C4) 



and u is the unpaired number density. 

Eqs. I|C3|I represent a formal solution only since w± 
must be determined from Eq. (|19[) in a self-consistent 
manner. This can be done by by combining Eqs. i|19|l . 
(|C1|) . and i|C3[l . It turns out that w+ and W- satisfy the 
same equation. Adopting the notation w± = w, we have 

2 UJqUJi[w 2 (uJq + UJl — 2uj q ) — 2wUJqUJi + 2ujQUJiUJq\ 2 

w 2 [w(ujq + uji) — 2uJoUJi] 2 (uJo — 2uj q )(uji — 2uj q ) 

(C5) 

with 

1 

and for the left side of Eq. 1C5|) 
V = V± = 



with 



5+ = 1- 



m = 



P±(T) 2 
I 

tanh 2 (e±/2fc B T) 
1 

w 2 (uj - UJl) 2 



(C6) 



2ujquji(uj — ujq)(uj — UJl) 
Eqs. I|C3|I now become 

\{uji-w)yi-x 2 -{2uj q A q y 

\/uj 2 (w - uji) 2 +uj 2 {w - UJ ) 2 



A d 

A^r 

A, 
A' 



\(w - uJp)U/i - x 2 - (2uj q A q ) 2 
\/uj 2 (w - uji) 2 + uj 2 {w - w ) 2 
(w - u) )(w - UJi) 



W 2 (uj — UJ q ) — WLOqUJ + UJqUjQUJl 
(w 2 ) — Auj q W + (Auj q UJ — UJqLl>i) 



W 2 (UJ — UJ q ) — WUJqUJ + UJ q UjQUJl 



(C7) 
(C8) 

(C9a) 
(C9b) 
(C9c) 
(C9d) 
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By introducing the variables 



where 



qo = 1 



2u q 



?1 = 1 



2m q 



and 



2 _ X — Gq _ go 
Xq ~ X-G, ~ q x 

Eq. I|C5|I can be greatly simplified: 



y qw — Qw + uj 



qow(w — Q) 



(CIO) 



(Cll) 



The signs "±" are for hole doping (a;, y > 0) and electron 
doping (x, y < 0), respectively. Obviously, this means 
that the solution has particle-hole symmetry: w depends 
only on the absolute value of y (or x) . For this reason, we 
shall regard y as the absolute value \y\ and ignore the ± 
sign in the following discussions. Solving the quadratic 
equation (|C11|) gives 



2 

<v) = (x-(x-G y 



(C12) 



- V (X - Gi)^Jx q (x q - y)(l - x q y), (C13) 

where r) = ±1. Eqs. (|CT2|) and l(€T3|) dete rmine the value 
of w and thus of w±. Since w + > in Eq. (|(J13[) . 
77 = +1 for w+ and 77 = —1 for u>_. 

Inserting Eq. (|C13|) into Eqs. (|U8ll and (|U9)l gives 



r(y) 



(Go — G\) 



2[G -e(y)][G 1 -e(y)] 



A q = ^nq(y) K = ^fnT(y) 



A,, 

with 

q(y) = - 



Go 



QS(y) X' 



i{y) 



fix 



[e(y) - Go][e(y) - d] x 
(y)2-[2e(y)-G -G 1 }x-G G 1 



(C14) 



(C15) 



(C16a) 



e(2/J - [2e(y) - G - Gijx - G Gi 

Eqs. I|C16|) can be further simplified by using Eq. HC10|) 
and recognizing that 

e(y)-G = -(x-Gi) ( V ab - a 2 ) 
e(y)~G 1 = (x-Gx) (b 2 - V ab) 
t(v)-X = -(X-Gi) (x q y + r)ab) 



= \ X q( X q~y) b = V 1 ~ X lV- 



(C17) 



In terms of a and b, Eqs. (|C14fl and l|C16(l can be rewrit- 
ten as 

r(y) = 



2a6 
a&x 

9(2/ = 

x q y 



T(y) = a 
S(y) = b . 





— a; 2 ) x 2 - 


-a 2 6 2 




(a 2 + 6 2 ) 


.X 2 




— x 2 ) x 2 - 


-a 2 6 2 



(a 2 + 6 2 ) x 2 



7(.) = (x-^) a2b2 " ( ^ )2 +X. 

Using Eq. (|C17() to convert a and b back to y yields 

m = (1 - 2/ 2 ) + (x g - y? 

2y/x q (x q -y)(l- x q y) 

q(y) = \J K - y)^? 1 - 2/) ^ 



T(y) = ^Jx(x q -y) g(y) 

s(y) = X/xixq 1 - 2/) 5(y) 

7(y) = (x-G 1 ) ^ (1 '^ ) +G 1 



9(2/) 



' x x q (Iy 2 — x 2 ) 
JJ xy 2 {l - 2x q y + x 2 ) ' 



Inserting the above results into Eqs. I|C15|) and taking 
into account the definition of x q in Eq. (|C10(1 gives Eqs. 
(El in Sect. V.A. 

A q = ^TrXjiXq 1 -y)( X q-y) 7" 

z y 



A d = -Tj-yjxixq 1 -y)g{y) 
GjCl 

Y 

_ (Gi-x)n fx 

A Xq I XqX 

2 \y 



A ^r = -^-\fx{xq - y)g(y) 



GxQ 



For T = 0, from Eq. (|C6|) . one obtains immediately 

y = y±= x 

because 5± — > 0. We thus have the results of Eqs. i|24|) 
discussed in Sect. IV. A. For T > 0, however, the self- 
consistency condition y = y + = y_ must be satisfied, 
which requires that 

y/l[l + S + T(y)} = - S.T(y)} 
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and thus that 6+ + 5- — 0. According to Eq. (|C7|) . this 
can be fulfilled only if 

_ tanh 2 (e + /2fc B T) +tanh 2 (e_/2yfc B T) 
2 ' 

which implies the existence of unpaired particles. The 
density of unpaired particles can be deduced from Eq. 
(El: 



This completes the derivation. Methods to obtain the 
actual gap values for given x and T from the above equa- 
tions have been discussed extensively in Section V. 



\2k B Tj \2k B Tj 
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